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NOMENCLATURE 


A 

a 

B(e) 

Bi 

Bi 

C(y) 

Di 

e 

Fi 

g(n,F^2) 

h 

h 

I 4 

kp 

M 

Ni 

s 

t 

'*^0 

^3j 

Vopt 

W(x) 


nuclear mass number, dimensionless 
oscillator parameter, fm 

average slope parameter of nucleon-nucleon scattering 
amplitude, fm^ 

defined in equation (35) 

defined in equation (10) 

correlation function, dimensionless 

defined in equation (11) 

two-nucleon kinetic energy in their center of mass frame, 
GeV 

defined in equation (12) 
defined in equation (25) 
defined in equation (41) 
defined in equation (42) 
defined in equation (43) 
defined in equation (44) 

Fermi momentum wave number, fm"^ 
nucleon mass, 938 MeV/c^ 
defined in equation (9) 
proton rms charge radius, fm 
defined in equation (13) 

average two-nucleon transition amplitude, MeV 
defined in equation (6) 

two-nucleon transition operator for nucleons 3 and j, MeV 
optical potential operator, MeV 
optical potential, MeV 
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X 


relative position vector of projectile, fm 
two-nucleon relative position vector, fm 
nuclear proton number, dimensionless 


y 

Z 


— > 

z 

collection of constituent relative coordinates for target, 
fm 

a(e) 

average ratio of real part to imaginary part of 
nucleon-nucleon scattering amplitude, dimensionless 

6 

defined in equation (23) 

6(x) 

Dirac delta function, dimensionless 

r 

defined in equation (47) 

Tl 

defined in equation (62) 

Y 

harmonic well charge density parameter, dimensionless 

T 

defined in equation (50) 

p 

nuclear density, fm"^ 

a(e) 

average nucleon-nucleon total cross section, fm^ or mb 

n j 

/ 

defined in equation (22) 

Slip 

defined in equation (48) 

o)T 

defined in equation (34) 

Subscripts; 


C 

charge 

m 

matter 

P 

projectile 

P 

proton 

T 

target 


nc Pauli correlations excluded 

pc Pauli correlations included 


Arrows over symbols indicate vectors. 
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ANALYTIC OPTICAL POTENTIALS FOR NUCLEON-NUCLEUS AND NUCLEUS-NUCLEUS 
COLLISIONS INVOLVING LIGHT AND MEDIUM NUCLEI 


Hari B. Bidasaria 
Old Dominion University 
and 

Lawrence W. Townsend 
NASA Langley Research Center 


SUMMARY 

Utilizing an optical model potential approximation to the exact 
nucleus-nucleus multiple-scattering series, optical potentials for nucleon- 
nucleus and nucleus-nucleus collisions are analytically derived. These 
expressions are applicable to light and medium cosmic ray nuclei as their 
single-particle density distributions are analytically determined, without 
approximation, from their actual harmonic well charge density distribu- 
tions. Pauli correlation effects are included through the use of a simple 
gaussian function to replace the usual expression obtained in the infinite 
nuclear matter approximation. 



INTRODUCTION 


Microscopic analyses of optical potentials for use in nucleon-nucleus 
and nucleus-nucleus scatterings typically utilize Woods-Saxon parameter iza- 
tions for the single-particle densities of the colliding systems. Although 
accurate as representations for the density distributions of most nuclei 
heavier than oxygen (ref. 1), they are only approximately correct for light 
and medium nuclei (3j<Z£8). For these nuclei, a more exact representation 
for the actual nuclear density is the harmonic well distribution (ref. 1) 
based upon the nuclear shell model. In addition to being a more accurate 
density representation, the harmonic well distribution, unlike the Woods- 
Saxon one, allows a completely analytic expression for the optical 
potential to be determined. In this paper, optical potentials for nucleon- 
nucleus and nucleus-nucleus scattering, with and without Pauli correla- 
tions, will be analytically derived for light and medium cosmic ray nuclei. 
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NUCLEON-NUCLEUS POTENTIALS 


For scattering between composite nuclear systems, the effective opti- 
cal model potential operator (ref. 2) is 

''opt = I t , , (1) 

6j « 

where tgj is the transition operator for scattering between the a 
constituent of the target and the j constituent of the projectile. From 
references 3-5, equation (1) gives a generalized expression for the optical 
potent i al 


W (x) = ApA.p /d^z p.|.{z) / d^y Pp(x + y + z) t(e,y) , (2) 

where t is the two-body transition amplitude averaged over the constitu- 
ent types (ref. 3) and p, and pp are the target and projectile 
single-particle number (matter) densities. 
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NUCLEON-NUCLEUS POTENTIALS 


Optical Potential Excluding Pauli Correlations 

Since equation (2) was derived using simple unsymmetrized product 
wavefunctions (refs. 3, 5), it does not include Pauli correlation effects. 
If the incident projectile is elementary (i.e., a nucleon), then Ap = 1 
and the projectile single-particle density is 

Pp(x + y + ?) = 6(? + ^ + 1) , (3) 

so that equation (2) reduces to 

W (x) = A.J. / d^z p^(z) t(e, X + z) , (4) 

where 

t(e, y) = -t^ exp[-y ^/2B(e)] , (5) 


and 


1 -3 

= (e/m) „(e) [c(e) + 1] [2. B(e)] . (6) 

In equation (6), e is the two-nucleon kinetic energy in their center of 
mass frame, and a(e), a(e), and 3(e) are the usual nucleon-nucleon scat- 
tering parameters (refs. 3-5). 

For light and medium target nuclei, the matter density is 


pj(r) = N.J. (C.p + Dj r^) exp(-F.|. r^) , (7) 

which results when the gaussian proton charge density is unfolded from the 
harmonic well charge distribution 

Pj(f) = Po [1 + (-f exp(-rW) , (8) 

using the general procedures given in references 4 and 5. In equation (7), 
the various parameters, written in terms of the charge density parameters, 
Y and a, are 


= p^a^/Ss^ , 

(9) 

= 1 + (3y/ 2) - (3Ya^/8s^) , 

(10) 
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( 11 ) 


7 d 

Dj = Ya^/ 16 s , 

Fj = l/4s^ , (12) 

where 

= (a^/4) - (rp^/6) , (13) 

and Tp « 0.87 fm (ref. 6) is the proton rms charge radius. Values for y 
and a are tabulated in ref. 1. The normalization constant, pg, is 
determined by requiring that 

4tt /” r^ p^(r) dr = 1 . (14) 

Incorporating equations (5) and (7) into equation (4) yields 

W(x) = t^ J d^z (C^ + z^) exp(-F.^ z^) ^ 5 ^ 

X exp[-(x + z)^/2B(e)] , 
which can be rewritten as 

W(x) = -2tt A.j.Nj t^ exp [-x^/2B(e)] /“ z^dz (Cj + Dj z^^ 

X exp[-(Fj + ig) z^] /q exp[-xz cos 0/B(e)] sinede . 

The angular integration is of the form 

Q exp(-qcos0) sin0d0 = [exp(q) - exp(-q)]/q , (17) 

giving 

W(x) = -2ir AjNj t^ exp[-x^/2B(e)] [B(e)/x] (18) 

X /” zdz (Cj + D.J. z^) exp{-[F.j. + (l/2B)]z^> 

X (exp[xz/B(e)] - exp[-xz/B(e)]} . 

The z integration can be rewritten as 

/“ dz (CjZ + Dj z^) exp{-[Fy + l/2B(e)]z^} {exp[xz/B(e)] (19) 
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-exp[-xz/B(e)]} = /“^ dz (C^z + D^z^) 

X exp{-[Fj + l/2B(e)]z^} exp[xz/B(e)] , 
which is of the general form (ref. 7) 

r / exp[-pz^ + 2qz] di = [(2l-")/pX’/p)^''^ (20) 


, d'’'^ [q exp(q^/p)] . 

Thus, the nucleon-nucleus optical potential, without Pauli correlations, is 
W(x) = (e/m)^/^ a(e) [a(e) + i] (21) 

X [Cj + (3fi-pD.|-/2) + x^] , 

where 

= (3 + , (22) 


and 


B = l/2B(e) . 


(23) 



NUCLEON-NUCLEUS POTENTIALS 


Optical Potential Including Pauli Correlations 
From ref. 8, Pauli correlation effects can be included by substituting 

Pil^i) Pi(r2) Pi(^l) Pi(^2) 

where 

g(rp r^) ” 1 " ^2p * 

with 

C(y) = 0.238 exp(-kp y^/5) . (26) 

In equation (25), kp = 1.36 fm“^ is the Fermi momentum wave number for 
spatially infinite ^nuclear matter. Incorporating the ansatz, equa- 
tion (24), into equation (2) yields a generalized optical potential, 
including Pauli correlations, of 

W(x) = ApA.p/d^z pj{z) /d^y Pp(x + y + z) t(e,y) [1 - C(y)] . (27) 

For an elementary projectile, this reduces to the nucleon-nucleus potential 
W(x) = A^ /d^z py(z) t(e, X + z) [1 - C( | x + z|)] , (28) 


which is separated into parts as 

W(x) = W^^(x) + Wp^(x) . (29) 

In equation (29), 


W^^(x) = A.|./d^z p.p(z) t(e, X + z) , (30) 

which was the potential in equation (4) for no correlations (nc). Its 
solution is given by equation (21). The second part is 

Wpc(x) = A.j./d^z p.|.(z) t(e, x + z) C(|x + z|) , (31) 

which is the Pauli correlation (pc) correction to the nucleon-nucleus 
optical potential. 

Incorporating equations (5), (7), and (26) into equation (31) gives 
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(32) 


Wp^(x) = 0.476 TT exp{-[(k^/5) + (l/2B(e))] x^} 

X /” z^dz (C.p + D.pZ^) exp{-[F.p + (l/2B(e)) 

+ (kp/5)] z^} /J exp{-[xz/B(e)) + (2kp xz/5)] cose} 

X sinede . 

The solution to equation (32) is obtained in the same manner as the 
solution to (21) since it is of the same general form. Hence, the Pauli 
correction to the potential is 

Wp^(x) " 0.238 (e/m)^^^ »{e) [c(e) + 1] (S (33) 

X [Cj + (3 to.j.Dy/2) + B? iDj DyX^] exp(-B. WjFjX^) , 

where 

= (Fj + B^)’^ , (34) 

and 

B^ = e + (kp/5) . (35) 

The nucleon-nucleus optical potential, with Pauli correlations, is 

W(x) = (e/m)^/^ a(e) [a(e) + i] 6^^^ [C^ (36) 

+ (3 fJjD.p/2) + D.|.e^fi^x^] exp(-eF.pn.|.x^) 

- 0.238 [C.p + (3 oj.pD.p/2) + B^ '^T 

2 

X exp(-B^o)^F^ X )} . 
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NUCLEUS-NUCLEUS OPTICAL POTENTIALS 


Optical Potential Excluding Pauli Correlations 

From equation (2), the generalized expression for the nucleus-nucleus 
optical potential is 

W(x) = ApA^ /d^z p^(z) fd^y Pp(x + y + z) t(e, y) , (37) 

where, as before, t is given by equations (5) and ( 6 ). For light and 
medium projectile nuclei, the matter density is 

Pp(r) = Np (Cp + Dp r^) exp(-Fp r^) , (38) 

where the parameters, Np, Cp, Dp and Fp are also given by 
equations (9) through (12). The target nuclei matter density is given in 
equation (7). 

With these assumptions, the optical potential is 

W(x) = -ApA^N.pNpt^ fdh (C.^ + D.^z^) exp(-FjZ^) (39) 


X fd^y [Cp + Dp (x + y + z)^] exp[-Fp (x + y + z)^] 

X exp[-y^/2B(e)] . 

Expanding the sums yields 

W (x) = -ApA^NpN^t^ {C.pCpI^ + D^CpI^ + C^Dpl 3 + D^DpI^} , (40) 


where 

I^ = /d^z exp(-F.pZ^) /d^y exp[-Fp (x + y + z)^] (41) 

X exp(-y^/2B(e)] , 

I^ = fz^d^z exp(-F.|.z^) fd^y exp[-Fp (x + y + z)^] (42) 

X exp[-y^/2B(e)] , 

I 3 = fd^z exp(-F.j.z^) fd^y (x + y + z)^ (43) 
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X exp[-Fp (x + y + z)^] exp[-y^/2B(e)] , 

I4 = /z^d^z exp(-F^z^) /d^y (x + y + z)^ (44) 

X exp[-Fp (x + y + z)^] exp[-y^/2B(e)] . 

Evaluation of Ii and I2 is straightforward since the angular depen- 
dencies are given by equation (17) and the spatial integrations are of the 
general form in equation (20). Thus, merely quoting the results, we have 

exp(-3FFpF^x^) , (45) 

and 

I 2 = TT^ (1.5 + B^Fp^r x^) exp(-grFpFj x^) , (46) 

where Oj is defined in equation (22) and 

r = (3Fp + 3F^ + F^Fp)"^ , (47) 

ilp = (3 + Fp)'^ . (48) 

The integral over y in I3 and I4 can be written as 

/d^y (x + y + z)^ exp[-Fp (x + y + z)^] exp(-gy^) (49) 

= /d^y (?^ + 2?ycos0 + y^) exp[-Fp (I + y)^] 

X exp(-ey^) , 

where 

I = X + z . (50) 

Two of the the three angular integrations in equation (49) are of the type 
whose solution is given by equation (17). The remaining angular integra- 
tion is of the form 

|g-q COS0 sin0d0 = -q"^ [exp(q) + exp(-q)] (51) 

0 
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+ q'^ [exp(q) - exp(-q)] . 

Hence, equation (49) becomes 

fd^y (x + y + z)^ exp[-Fp (x + y + z)^] exp(-gy^) (52) 

= [1.5 + fip$^ (x + z)^] exp[-gfipFp (x + z)^ . 

Inserting equation (52) into I 3 and I 4 and using equations (17), (20), 
and (51) to evaluate the resulting integral forms gives 

1 3 = TT^ (1.5 r x^) exp(- 6 F^Fpr x^) , (53) 

and 

1 4 = {2.25 + 3 . 750 ^r + [1.5f?^r (Fp + F^) (54) 

- 5e^FpF-pr^] x^ + e^F^F^r^ x"^} exp(-BF-pFpr x^) . 

Upon combining these integrals, the nucleus-nucleus optical potential, 
without Pauli correlations, is 

W(x) = -ApA-pNpN^ (e/m)^^^ o(e) [a(e) + i] (ngr)^/^ (55) 

X ^C^Cp + D-j-Cpr (1.5 Jlp^ + e^FpT x^) + C^DpT 
X (1.5 + e^F^r x^) + DpD.j.r { 2.25 + 3.75B^r 

+ [1.56^r (Fp + Fj) - 53^FpF.j.r^] x^ + B^F^F^r^ x^} 

p 

X exp(-FpF.|-rB X ) . 
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NUCLEUS-NUCLEUS OPTICAL POTENTIALS 


Optical Potential Including Pauli Correlations 

From equation (27), the nucleus-nucleus optical potential, including 
Pauli correlations, is 


3 "^ 3 "^ 

W(x) = A A /d z p (z) /d y p (x + y + z) t(e, y) (56) 

P T T P 

X [1 - C(y)] , 

with C(y) given ^by equation (27), p-j- by equation (7), pp by 

equation (38), and T(e,y) by equation (5). Separating the potential as 
in equation (29), we find 

W(x) = W^^(x) + Wp^(x) , (57) 

where Wpq(x) is 

Wnc(x) = ApAj Jd^z p-|.(z) /d^y Pp(x + y + z) t(e, y) , (58) 

whose solution is given in equation (55), and W (x) is 

pc 

Wp^(x) = -ApA^ /d^z p^ (z) /d^y Pp (x + y + z) (59) 


X t(e, y) C(y) . 

Inserting equations (5), (7), (27), and (38) into equation (59) yields 
Wp^(x) = 0.238 ApA^NpN^t^ /d^z (C.^ + D^^z^) exp(-F^z^) (60) 


X /d^y [Cp + Dp (x + y + z)^] exp [-Fp (x + y + z)^] 

X exp(-ey^) exp(-kpy^/5) . 

This expression can be put into the same form as equation (39) by using 
equation (35) to combine the exponentials in y. Hence, the solution, 
which is of the form given by equation (55), is 



X C^Cp + D^CpF^ (1.5 0)p^ + x^) 

X C^DpF^ (1.5 uij^ + X^) + D^DpF^ 

X {2.25 + 3.75 + [1.5 (Fp + F^) 

- 5B^3 ^2 , g^4p2p2,3 ^4j 

X Gxp ( - Fp FjF^ B^^ X ) 

where 


= (Bj^Fp + B^Fj + FjFp) 


-1 


(62) 


Combining the results in equations (55) and (61) gives a nucleus-nucleus 
optical potential, including Pauli correlations, of 


W(x) = -ApAyNpNj (e/m)^^^ a(e) [a(e) + i] 


^znp/z 


(63) 


X ^C.|.Cp + DyCpF (1.5 fip"^ + B^Fp^r x^) + C.j.Dpr 

X (1.5n^^ + B^F^^r x^) 


+ DpD^r {2.25 + 3.75B^r + [1.53^r(Fp + F.^) - 5 B^FpF^r^]x^ 

+ 3^FpF-^^r^ x^}^ exp(-FpF.pr3x^) - 0.238 
X ^C^Cp + D^CpF^ (1.5 0)-^ + B^^Fp^r^x^) 

+ CyDpF^ (1.5 (o“^ + B^^Fy^r^x^) + DpD.j.r^ 

X {2.25 + 3.75 B^^r^ + [1.5 B^^r^ (Fp + F.^) - 5B^^FpF^Fj^^]x^ 

+ B^^^Fp^F/r^^ x"^}) exp (-FpF^-r^B^x^) . 


13 



EXAMPLES 


To illustrate these results, the real and imaginary parts of the 
optical potentials for proton-carbon and carbon-carbon scattering are 
plotted in figures 1 through 4. The incident projectile kinetic energy is 
200 MeV/nucleon. Values for the carbon charge density parameters, a and 
Y, were taken from reference 1. Values for the nucleon-nucleon scattering 
parameters, a(e), B{e), and a(e) were taken from the compilations in ref- 
erences 9 and 10. 

In figure 1 real parts of the proton-carbon optical potential are 
plotted, with and without Pauli correlations, as determined from 
equations (21) and (36). Note that the Pauli correlations reduce both the 
depth (strength) and width (range) of the interaction. Similar results are 
shown in figure 2 for the imaginary parts of the proton-carbon potential. 
The strength (depth) of the imaginary part of the potential is also greater 
than that for the real part. 

Figures 3 and 4 display comparable results for carbon-carbon scatter- 
ing. Again, the increased strength of the imaginary part of the interac- 
tion to the real part is evident, as is the reduced depth and width when 
Pauli effects are included. Note also that the well depth is much greater 
than in the proton-carbon interaction. This is indicative of the increased 
number of possible interactions in nucleus-nucleus collisions. 
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CONCLUDING REMARKS 


In summary, analytic expressions for nucleon-nucleus and nucleus- 
nucleus optical potentials, applicable to light and medium cosmic ray 
nuclei, have been determined. Pauli correlation effects were included in 
an approximate, yet analytic, manner and resulted in a reduction of both 
the strength and range of the overall interaction. The use of these 
expressions, at energies below 100-200 MeV/nucleon, should be carefully 
monitored since the increased importance of spin-dependent effects (which 
are ignored here) with decreasing energy, may render these expressions only 
crudely approximate at best. Above 400 MeV/nucleon, the results herein may 
be considered quite accurate for most calculations. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 1982 
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Figure 1. The real part of the proton -carbon 
optical potential at 200 MeV incident proton 
kinetic energy. 
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Figure 2. Imaginary part of the proton -carbon 
optical potential at 200 MeV incident proton 
kinetic energy. 
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Figure 3. Real part of the carbon-carbon 
optical potential at 200 MeV/ nucleon in- 
cident kinetic energy. 
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Figure 4. Imaginary part of the carbon-carbon 
optical potential at 200 MeV/ nucleon incident 
kinetic energy. 
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